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Abstract. We consider an open system of non-interacting electrons consisting of a 
small sample connected to several reservoirs and temperature or voltage probes. We 
study the non-linear system of equations that determines the probe parameters. We 
show that it has a unique solution, which can be computed with a fast converging 
iterative algorithm. We illustrate our method with two well-known models: the 
three-terminal system and the open Aharovov-Bohm interferometer. 



1 Introduction 

Thermodynamic quantities such as entropy, temperature, or chemical potential play a funda- 
mental role in our understanding of equilibrium phenomena. They are given sound microscopic 
meanings within the framework of equilibrium statistical mechanics. The concept of local ther- 
mal equilibrium allows us in principle to define these quantities in interacting systems close to 
equilibrium. However, extending these definitions far from equilibrium and/or to non-interacting 
systems where local equilibrium does not make sense is a much more delicate issue (see the 
discussions in [1, 2, 3]). In this paper, we shall consider an operational point of view, giving to 
local intensive parameters the values measured by external probes. 
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Such an experimental approach is well-known in the mesoscopic community: in the description 
of electric transport in a multi-terminal system, in which all the terminals have the same tempera- 
ture (typically T = 0), one often introduces a voltage probe [4] to sense the local electrochemical 
potential by connecting an additional electronic reservoir under the zero electric current condi- 
tion: the chemical potential of the probe is tuned so that there is no net average electric current 
into it. In the same spirit, setting all the terminals to the same chemical potential, a temperature 
probe is obtained by requiring that the temperature of the corresponding reservoir is tuned such 
that there is no average heat current into it [5]. 

In the scattering approach of Landauer and Biittiker (see Section 2), the existence and uniqueness 
of such parameters are usually accepted on physical grounds, but we think it is important and 
interesting to obtain a rigorous mathematical foundation for these fundamental parameters. In the 
linear response regime, a rigorous proof has recently been given [6]. Here, we shall extend these 
results to the far- from-equilibrium regime and furthermore provide an efficient numerical method 
for computing their values. More explicitly, this paper is organized as follows: in Section 2, we 
describe the framework, in Sections 3 and 4, we present our main results and their proofs. Finally, 
in Section 5, we illustrate our method by considering two well-known models: the three-terminal 
system and the open Aharovov-Bohm interferometer. 
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2 Framework 

We consider a multi-terminal mesoscopic system, that is, a small system S connected through 
leads to several infinitely extended particle reservoirs (see Figure 1). We assume that the trans- 
port properties of this system can be described within the Landauer-Biittiker framework. More 
precisely, we consider iV reservoirs in equilibrium at inverse temperature = l/^B^i) and 
chemical potential = 1, . . . , N). The corresponding Fermi-Dirac distributions are 

fi(E) = f(E, A, IH) = [l + e^ E -^}-\ (1) 

For simplicity, here and in what follows we set the Boltzmann constant k&, the Planck constant 
h and the elementary charge e to unity. 

In the Landauer-Biittiker formalism, one neglects all interactions among the particles and con- 
siders the small system as a scatterer for the particles emitted by the reservoirs. Thus, the 
small system is completely characterized by the one-particle on-shell scattering matrix S(E) = 
[Sij ; mn(E)], where the indices i, j 6 {1, . . . , N} label the outgoing/incoming terminals and for 
each pair the indices m E {1, . . . , Mi(E)} and n E {1, . . . , Mj(E)} label the open chan- 
nels in terminals i and j, respectively. The matrix element Sij- mn (E) is the probability amplitude 
for a particle with energy E incident in channel n of terminal j to be transmitted into channel m 
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of terminal i. The corresponding total transmission probability Uj(E) that a particle with energy 
E goes from terminal j to terminal i is given by [7] 

Mi(E) Mj(E) 

tiAE)= £ £ |% ;m n(i?)| 2 . (2) 
m=l n=l 

The unitarity of the scattering matrix immediately yields the following identities, 

N N 

£%(£) = Mj(E), £%(£) =M i (£). (3) 
i=i j=i 




Figure 1: A multi-terminal system: The sample is connected through leads to N = N p + N t 
particle reservoirs. The reservoirs 1, . . . ,N p are probes and the reservoirs N p + 1, . . . , N p + N t 
are thermostats driving the system out of equilibrium. 

The expected stationary electric and heat currents in lead i e {1, . . . , iV} are given by the cele- 
brated Landauer-Biittiker formulas [8, 9], 

tij{E)fj{E)] dE, (4) 

Uj{E)fj{E)] (E — fMi)dE. (5) 

From a phenomenological point of view, these expressions can be easily understood: tji(E)fi(E) 
is the average number of particles with energy E that are transmitted from terminal i to terminal 
j, and tij(E)fj(E) is the same but from terminal j to terminal i. Therefore, Jj (Jj) is the net 
average electric (heat) current in lead i, counted positively from the i-th terminal to the system. 
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Mathematical derivations of these formulas (including existence of a stationary regime) rest on 
the assumption that the leads are infinitely extended and act as reservoirs [10, 11]. 

Considerable interest has been devoted to electric transport in which all the terminals have the 
same temperature. In this context, an important concept emerged: the voltage probe [4]. A 
voltage probe is a large physical component used in mesoscopic experiments to sense the local 
electrochemical potential. Theoretically, such a probe is modeled as a reservoir under the zero 
electric current condition: the chemical potential of the probe is tuned so that there is no net 
average electric current into it. If all the terminals have the same temperature, then in general 
there will also be a heat current into the probe. In this case, we will consider this heat current 
as dissipation. In the same spirit, setting all the terminals to the same chemical potential, a 
temperature probe is obtained by requiring that the temperature of the corresponding reservoir is 
tuned such that there is no average heat current into it. Note that in this case, there may be some 
charge dissipation into the temperature probe. 

Let us decompose the N terminals as follows: the first N p reservoirs are temperature or volt- 
age probes and the remaining N t reservoirs are the thermostats maintaining the system out of 
equilibrium (see Figure 1). In the voltage probe configuration, all the reservoirs are at the 
same inverse temperature (3 = (3\ = ■ ■ ■ = (3n, the chemical potentials of the thermostats 
Pt = (l^Np+i, ■ ■ ■ >/A/Vp+jVt) are given and we have to determine the probe parameters /7 P = 
(pi, . . . , /ijv p ) such that I p = (Ji, ... , Jjv ) = 0. Similarly, in the temperature probe configura- 
tion, all reservoirs are at the same chemical potential p = \i\ = ■ ■ ■ = p N , the thermostat inverse 
temperatures f3 t = (Pn p +i, • • • , PN P +N t ) are fixed and we have to determine the probe parameters 
P P = {fix, . . . , /3 Np ) in order to satisfy J p = ( J u . . . , J Np ) = 0. 

To our knowledge, no result is available on these two problems beyond the linear approximation 
around global equilibrium (i.e., linear response theory, see [5]). The same remark applies to other 
approaches to the determination of local intensive thermodynamic parameters (see e.g., [12, 13]). 

3 Results 

Note that in both configurations the self-consistency condition 

J p (/3, p t ; ftp) = 0, or J p (p, /3 t ; f3 p ) = 0, 

is a system of iV p non-linear equations with iVp unknown. From a mathematical perspective, it 
is not at all obvious that such a system admits a solution. Moreover, if a solution exists, it may 
not be unique. Our main result ensures existence and uniqueness of reasonable solutions to these 
equations. 

We shall make the following general assumptions on the lead Hamiltonians and scattering matrix: 

(A) There exists a constant E such that Mj(E) = for all E < E and j e {1,...,N}. 

(B) Mj(E) < C(l + \E\)i for some constants C and rj and all E and j G {!,..., N}. 
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(C) For every j E {1, . . . ,N P } there exists a set £j C E of positive Lebesgue measure such that 

N p +Nt 

M 3 {E)= t ij (E)>0, 

i=N p +l 

for all E G £j. 

Condition (A) merely asserts that the lead Hamiltonians are bounded below. Condition (B) limit 
the growth of the number of open scattering channels as function of the energy and is satisfied 
by any physically reasonable lead Hamiltonian. Finally, Condition (C) can be roughly rephrased 
as follows: any probe is connected through an open scattering channel to some thermostat. 

To formulate our main result, let us denote 

/i = min{/i7v p+ i, . . . , /ijv p+ 7v t }, 
\i = max{fi Np+ i, . . . , fiN p +Nt}, 

the minimal/maximal chemical potential of the thermostats and define in the same way (3 and (3. 
Theorem 1. Under the above assumptions, the following hold: 

(1) The self-consistency condition I p ((3, (2 t ; /%,) = has a unique solution jl p = /2 p (/3, £t t ) m the 
set {(//!,..., fx Np ) \nj E [//,//]}. 

(2) The self-consistency condition J p (/i, (3 t ; (3 P ) = has a unique solution (3 P = (3 p (fi, (3 t ) in the 

set{(p l ,...,p Nv )\p j e\p,M 

(3) In both cases the solution can be computed by means of a rapidly convergent algorithm (see 
the next sections for details). 

Remarks. 1. The restriction on the solution is physically reasonable. We do not expect a 
temperature probe to measure a value below the smallest thermostat temperature or above the 
highest one. The same remark applies to voltage probes. 

2. An alternative approach to probing local intensive parameters is to adjust both (3 P and /2 p in 
such a way that the electric and heat currents vanish: I p — J p — 0. Such probes thus measure 
simultaneously the temperature and the chemical potential. Note that in this case there is no 
dissipation at all into the probes. Our method does not apply directly to this situation, basically 
because the function f(E, (3, fi) — f(E, /?', //) does not preserve its sign as E varies if f3 ^ f3' 
and |i / To our knowledge, no result is available for such dual probes beyond the linear 
approximation around global equilibrium (see [14, 6]). 
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4 Proofs 

Let us discuss first the voltage probe configuration. Using the relations (3), the self-consistency 
condition may be written as 

- r 

h = J2 f( E ^^j)[ M A E ) S ij ~ tij(E)]dE = 0, (6) 
for i = 1, . . . , N p . Under Assumptions (A) and (B), 

H^X j (ji) = f f(E,P,fi)M 3 (E)dE, (7) 

defines a strictly increasing continuous function. We shall denote by X ^ Pj{X) the reciprocal 
function. The key idea of our approach is to work with the variable X = (X 1 (/i 1 ), . . . , X Np (p Np )) 
instead of /2 p . 

Let F : R Np -> R Np be defined as 



F(X) = J^ j f{E,0, N {Xj))Uj{E)dE 

N p +Nt 

+ / f{E,{3,H)UAE)&E. 

„■ AT i 1 ** 



j=N p +l ■ 

Then, we can rewrite the self-consistency condition (6) as a fixed point equation 

F(X) = X. (8) 
Set X_j = Xj(ii), Xj = Xj(jt) and denote 

E = {X = (X 1; . . -,X Np ) | Xj G [X^Xj]}. 
Notice that the condition X G E is equivalent to fij G [/a, ju] for all j G {1, . . . , X p }. 
Lemma 2. F(E) c S. 

Prao/ Let X G E. The monotony of p i-» f(E, (3, /i) implies f(E, (3, Pj{Xj)) < f(E, 0, JJ) for 
j = 1, . . . , X p and /(£?, /3, pj) < f(E, /3, /I) for j = X p + 1, . . . , X p + X t . The identities (3) 
yields 

Fi(X)< [ f(E,P,-p)Mi(E)dE = Xi. 



Proceeding similarly, one shows 

F(X)> I f(E,p,/i)M i (E)dE = X i 



□ 
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Under Conditions (A) and (B), the function F is continuous. Since E is compact and convex, it 
follows from Lemma 2 and the Brouwer fixed point theorem that (8) has a solutions X* G E. 
We shall use Condition (C) to ensure uniqueness of this solution. In the next lemma, we use the 

norm||X|| = J^fZi \Xj{ 

Lemma 3. Under Assumptions (A), (B) and (C) there exists a constant 9 < 1 such that 

\\F(X) -F(X')\\ < 6\\X-X'\\, (9) 

for any X,X' G E. 

Proof. Denote by D(X) the derivative of the map F at X. Then one has 

F(X)-F(X')= [ D(tX + (1 -t)X')(X -X')dt. 
Jo 

Since E is convex, the estimate (9) holds for any X,X' 6 E with 

9 = max\\D(X)\\, 

where the matrix norm is given by 

\\D(X)\\ = max_J2\ D iAX) 

A simple calculation yields 



i<j<N P *r , 



where the function g(E, (3, fj) = d^f(E, (3, /i) is strictly positive. It follows that 

/ Mj(E)g(E, (3, Hj(Xj))dE 



< i - 



i=i 

and hence 

9 < 1 — min 



/M,(^( J E,/3,/i,(X J ))dE' 
fMj(E)g(E,f3,Li)dE 



^ j M 3 (E)g(E^^)dE- 



Condition (C) clearly implies that 6 < 1 . □ 

It follows from Lemmas 2, 3 and the Banach fixed point theorem that (8) has a unique solution 
X* in E. Moreover, the sequence of iterates X n = F(X n _i) converges to X* for any initial 
value X G E with the estimate 

\\x n -x*\\ < e n \\x -x*\\. 
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In the temperature probe configuration, one may proceed in a completely similar way in terms of 
the functions 

P^Y j (0) = J f&frME-riMjWdE, (10) 
their reciprocal Y i-» /3j(Y) and 



i=i J 

N p +Nt 

„■ AT r 1 *f 



j=N p + l ■ 

A natural set £ can then be defined as before. The crucial observation is that the function 

dpf(E, f3, fi)(E — ii) has a constant sign. 

Remark. Strictly speaking, Lemma 3 does not hold at zero temperature because the Fermi 
function / is not positive in this case. Nevertheless, one easily shows that, under Assumptions 
(A)-(C), the estimate 

\\F{X)-F{X')\\<\\X-X% 

holds for 1,1' 6 S,I / X' provided [//,//] C C\j£j. The uniqueness of the fixed point X* 
immediately follows. Moreover, it also follows that the sequence of iterates X n = F(X n _i) 
converges to X* for any choice of X G S, although without a priori control on the speed of 
convergence. 



5 Examples 

As a first example, let us consider the one-channel three-terminal system represented in Figure 2, 
where two thermostats (2 and 3) drive the system (a perfect lead) out of equilibrium and a probe 
(1) is connected to the system by a 3 x 3 scattering matrix S. 

Let us consider the energy-independent scattering matrix introduced in [15]: 

f-(a + b) 

S=\ a b , (11) 

\ \fe b a ) 

where a = |(\A — % £ — 1), b = — 2e + 1) and e G (0, |]. Here, e = corresponds to 

the uncoupled situation (which is excluded) and e = | to the maximally coupled one. Let us set 
T = Ti = T 2 = T 3 and define the energy interval in (4)-(5) as [0, oo). If T = 0, then in the 
linear regime one can compute analytically the self-consistent parameter /i* [7, 14]: 
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Figure 2: A one-channel system with two thermostats (2 and 3) and one probe (1). 



We have checked that our numerical results are consistent with the relation (12). In the non- 
linear regime, we made the following observations: Let T > 0, //2,£*3 G IR be fixed, then the 
sequence {F n (X )}^ =0 , with X G E, converges and gives rise to a value nl independent of e 
and conveniently written as 



V*{T, /i2, M3, e) = ^ 2 ~t^ 3 + A/"(T, /i 2 , /i 3 ) 



(13) 



where the function M(T, /i 2 , /i 3 ) measures the non-linearity. Note, in particular, that the weak 
coupling limit e — > does not lead to a different value of n\. Since, by Theorem 1, 

jU*(T,//2,jU 3 ,e) G [min{/i 2 , /i 3 }, max{/i 2 , /x 3 }], 

one deduces that A/"(T, /i 2 , /i 3 ) 6 [— A/i/2, A/i/2], with A/x = |/i 2 — /i 3 |. In Figures 3 and 4, 
we have plotted the temperature and potential dependence of Af(T, /i 2 , /i 3 ), respectively. Let us 
recall that Af(T, /i 2 , /i 3 ) = corresponds to the linear case. 

Note that the curve in Figure 4 reaches a constant value TV^T) = lim^^oo J\f(T, /i 2 = 0, // 3 ) 
as /i 3 increases. Interestingly, this is also the case for other values of T and we observed the 
following scaling law: 

Afoo(XT) = Woo(T) , VAX). (14) 

If we attach more probes to the lead and describe all the connection points in terms of the same 
scattering matrix 5* (see [14] for the construction of the global scattering matrix), then we found 
that all the probes measure the same value, as if all the probes were connected to the same 
point, but in general this value does not coincide with the one-probe measurement (since adding 
more probes somehow perturbs the system). This phenomenon can be easily understood: for 
example, if two probes are attached to the lead, then one can compute analytically the global 
4x4 transmission matrix {Uj(E)}, and one finds that it is symmetric and that tsi(E) = t 32 (E) = 
t±i(E) = t 4 2(E). This means that the two probes are equally coupled to the left and right 
thermostats and, consequently, that ji[ = ji* 2 . Note, however, that this is not true in general. 
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Figure 3: The temperature dependence of Af(T, /x 2 = 0, /x 3 = 100). 




Figure 4: The potential dependence of M{T = 1, fi 2 = 0, 7/3). 



10 



Temperature and voltage probes far from equilibrium 



As a second example, we consider an Aharonov-Bohm (AB) ring threaded by a magnetic flux 
<f) and with a quantum dot (QD) embedded in one of its arms. This system has been subject to 
intensive investigations both in the independent electron approximation [12, 16] and including 
interaction effects [17, 19, 18]. We shall study a discrete (tight binding) independent electron 
model closely related to the work by Aharony et al. [20] (see Figure 5). This theoretical model is 
supposed to imitate an experimental setup [ ]. It is assumed that a gate voltage V is applied on 
the QD, allowing to vary the energies of its eigenstates. Let us write tq D = y^Tq^e 101 ^ for the 




Figure 5: The open AB interferometer. A magnetic flux <\> crosses the ring and a QD is placed in 
the upper branch of the ring. The terminals 1 to 12 are voltage probes and the terminals L and 
R are thermostats, hr, kp and ko are the hopping constants coupling the ring to the thermostats, 
the probes and the QD. The energy levels of the dot are V, V + U and V + 2U. 



transmission amplitude of the QD. At fixed energy E, the total transmission probability t LR from 
the reservoirs L to R depends on the gate voltage V and is a periodic function of the AB-flux (p. 
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Expanding this function as a Fourier series one gets: 



t LR (V, 4>) = A(V) + B(V) cos(0 + 0(V)) + ■■■ . 



(15) 



It is well known that in the absence of dissipation (i.e., for a closed interferometer in the terminol- 
ogy of [20]) the Onsager-Casimir reciprocity relations [ ] imply that the phase (3(V) can only 
take the values and it. Hence, as the gate voltage V varies, the phase f3(V) makes abrupt jumps 
between these two values. However, dissipation can change this picture. By adding purely ab- 
sorbing reservoirs (i.e., allowing only outgoing currents) along the branches of the ring Aharony 
et al. [20] found criteria as to when the "experimental" phase /3(V), which depends on the details 
of the opening (i.e., the coupling the absorbing reservoirs), is a good approximation of the "in- 
trinsic" phase ctQ D (V) of the QD. Here we present some numerical results showing that one may 
capture the main properties of o;qd without introducing any charge dissipation in the absorbing 
reservoirs, i.e., that (3 behaves essentially as «qd even if one replaces the absorbing reservoirs of 
[20] by voltage probes, which we recall allow only heat dissipation. However, instead of consid- 



Figure 6: The "intrinsic" phase «qd (blue dashed line) and the "experimental" phase (3 (red 
solid line) as a function of the gate voltage V applied on the QD. This corresponds to the setup 
of Figure 5 with a three-level QD, all temperatures equal to zero, the chemical potential of the 
thermostats are /jll = and fx R = 0.2. The couplings to the thermostats, probes and QD are 
kr = 0.5, kp = 0.5 kx, and kp, = 0.01 kp, respectively (the parameter kp plays a role similar to 
J x in [20]). All remaining parameters are set as in [20]. 

ering the expansion (15), we tried to be closer to actual experimental measurements by extracting 
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the "experimental" phase (3 from the Fourier expansion of the steady electric current between the 
two thermostats: 

I L = -I R = I (V) + h(V) cos(0 + /3(V)) + ■■■ . 

The results are shown on Figure 6. One sees that the curve (3(V) follows relatively closely 
«q D (V), and in particular reproduces accurately the successive jumps of «qd(V) from 1 to 
(the values have been normalized, thus 1 corresponds to n in the paper [20]). 
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